Interplay between disorder and local field effects in photonic crystal waveguides 
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We introduce a theory to describe disorder-induced scattering in photonic crystal waveguides, 
specifically addressing the influence of local field effects and scattering within high-index-contrast 
perturbations. Local field effects are shown to increase the predicted disorder-induced scattering 
loss and result in significant resonance shifts of the waveguide mode. We demonstrate that two types 
of frequency shifts can be expected, a mean frequency shift and a RMS frequency shift, both acting 
in concert to blueshift and broaden the nominal band structure. For a representative waveguide, we 
predict substantial meV frequency shifts and band structure broadening for a telecommunications 
operating frequency, even for state of the art fabrication. The disorder-induced broadening is found 
to increase as the propagation frequency approaches the slow light regime (mode edge) due to 
restructuring of the electric field distribution. These findings have a dramatic impact on high- 
index-contrast nanoscale waveguides, and, for photonic crystal waveguides, suggest that the nominal 
slow-light mode edge may not even exist. Furthermore, our results shed new light on why it has 
hitherto been impossible to observe the very slow light regime for photonic crystal waveguides. 

PACS numbers: 42.70.Qs, 42.25.Fx, 42.79.Gn, 41.20.Jb 
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Photonic crystal (PC) structures comprised of high- 
index-contrast cavities and waveguides offer a rich de- 
gree of control over light-matter interactions, leading to 
trapped 1 and slow light modes - ' buried within a forbid- 
den photonic band gap. In a planar PC semiconductor 
system, waveguide modes can be completely bound below 
the light line. However, manufacturing imperfections re- 
sult in fabrication disorder that breaks the translational 
invariance of a nominally perfect lattice, and causes ex- 
ternal scattering of the bound modes. From a theoretical 
perspective, the role and impact of fabrication disorder 
on the effect of slow-light slab waveguides is becoming 
better understood' -5 ; minute sidcwall imperfections act 
to perturb the propagating mode causing the light to out 
scatter through radiation modes or backscatter and mul- 
tiple scatter within the waveguide (, ~' s . Scattering losses 
become particularly pronounced in the slow light regime 
as the local density of states of the mode into which 
light can scatter increases. Although many experiments 
have confirmed this slow light loss behaviour , an open 
question that still remains is what is the effect of disorder 
on the band structure? The answer to this question in- 
volves a complex interplay between local field corrections 
and enforcing boundary conditions between parallel and 
perpendicular field components. Experimental evidence 
for a dramatic reduction of transmission near the mode 
edge is observed to take place much sooner than that pre- 
dicted for the scaling of scattering loss, e.g., see Ref. 6. 
This suggests that some unknown effect is either shift- 
ing the mode edge or broadening the band structure (see 
Fig. 1(a)). Here we address this question directly and 
show that surprisingly large changes in the band struc- 
ture occur which have not been anticipated before. Our 
findings introduce a dramatic revision of our present un- 
derstanding of the role of disorder-induced scattering in 
these intriguing nanostructures. 

The usual theoretical approaches to modelling 
disorder-induced scattering, employ standard perturba- 



tion theory where the nominal (disorder free) electric 
field is used with the dielectric index change to model 
polarization scattercrs, through P^is = AeEo, with Eo 
the unperturbed field and Ae is the dielectric contrast 
change resulting from a spatial perturbation (e.g., air to 
silicon gives Ae « 12) 12 . For high- index-contrast pertur- 
bations, as shown by Johnson et al. , this polarization 
scatterer is problematic for several reasons: (i) the paral- 
lel components of the electric field and the perpendicular 
components of the displacement field must be continuous 
across the surface, and (ii) the index change results in lo- 
cal field corrections. Andreani and Gerace 14 attempted 
to estimate the magnitude of this error in their theory by 
comparing with a simple numerical supercell calculation 
of perfect hole shapes with different radii, and they con- 
cluded that the effects of local fields was not important. 
Wang et al. employed the well known slowly-varying 
surface approximation 1 \ valid only for smooth bumps, 
and demonstrated the impact on increasing the scattering 
losses. For PC cavity systems, Ramunno and Hughes" 1 
showed that a quickly-varying surface perturbation can 
cause mean resonance shifts of a strongly confined cav- 
ity mode 11 '. RMS ensemble average shifts and band edge 
broadening have also been predicted (without including 
local field effects) for intrinsically lossy coupled-cavity 
PC waveguides 1 '' 18 . 

With regards to PC waveguides, to the authors' knowl- 
edge, there have been no calculations nor any awareness 
of disorder-induced resonance shifts. Neither has anyone 
computed disorder-induced losses in the presence of real- 
istic rapidly- varying surface perturbation by properly ad- 
dressing the two problematic criteria above. In this work, 
we overcome this limitation, and introduce a straightfor- 
ward optical scattering theory that allows one to compute 
scattering losses and resonance shifts within an ensem- 
ble averaging procedure, while accounting for local field 
effects. The ensemble average statistical calculation ap- 
plies to measurements of an ensemble of nominally iden- 
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FIG. 1: (Color online) (a) Nominal dispersion of PC wave- 
guide mode (blue), and the broadened disorder-induced band 
structure that we will introduce later (grey shading), (b) 
Electric field Bloch mode near the band edge along symmetry 
planes in the unit cell, (c) Example of a disordered hole profile 
that satisfies the statistics used in the calculation (blue). The 
ideal radius (dashed black) and correlation length (short red 
arc) are shown for reference, where R indicates the nominal 
radius of the unperturbed hole and a is the RMS fluctuation. 



tical structures, and (or) it applies to the statistics of one 
waveguide which contains many nominally identical unit 
cells. This is a reasonable first model to report, since 
typical PC waveguides of only 1 mm-length contain al- 
ready many thousands of unit cells; our methodology can 
also help guide future coherent scattering theories, where 
such effects were neglected' . Using representative calcu- 
lations, we subsequently demonstrate that significant and 
unusual disorder-induced changes in the band structure 
can occur, as well as significant modifications (increases) 
to the predicted scattering losses. These predictions are 
interesting in their own right and are important for the 
analysis and interpretation of related experiments. 

The total electric field, in the presence of disorder, can 
be calculated from an integral solution of the Maxwell 
equations, through 

E(r;w) =Ei(r;w)+ f dr' G(r, r'; w) • P(r'; w), (1) 

J all space 

where P(r';w) is the polarization-like density due to the 
disorder in the system, Ej(r;u;) is the electric field in 
the ideal (no disorder) system 12 , and G(r, r';w) is the 
photon Green function. The Green function is simply 
a dipole solution to the Maxwell wave equation and it 
contains information about how light scatters as well as 
the local photon density of states (i.e., LDOS(r,uj) oc 
Trace{Im[G(r, r; w)]}). For convenience, we partition the 



Green function into contributions from the bound wave- 
guide mode, radiation modes above the light line, and 
other modes as G(r, r';ui) = Gb(i", r';cj) + GR(r, r';w) + 
Go(r,r';w). The bound mode Green function is given 
analytically from properties of the bound mode' 3 ' 19 : 
G B (r,r';w) = if^-[e fe (r) ® e* k (r r y k <- x - x '*>0{x - x 1 ) + 

e* k (r) (g) efe(r') e %k ^ x ~ x ^Q{x' - a;)], where v g is the group 
velocity, efc(r) is the Bloch mode electric field normal- 
ized by J cell dre(r) |efc(r)| 2 = 1, ® is a tensor product, 
e_fc(r) = et(r), and <d(x) is the Heaviside step func- 
tion. For our calculations below, we use a Wl PC wave- 
guide formed from a semiconductor membrane by omit- 
ting a row of holes in a two-dimensional triangular array 
of holes. The lattice pitch is a = 480 nm, the hole ra- 
dius r — 95 nm and the slab thickness is h = 160 nm. 
The dispersion of the ideal waveguide mode is shown in 
Fig. 1(a). The electric field Bloch mode distribution is 
shown on symmetry planes in the unit cell in Fig. 1(b) 
for a wave vector near the mode edge. 

The disorder fluctuations of interest must closely cor- 
respond to real images of fabricated PC waveguides 2 ", 
and thus we consider disorder in PC slab structures that 
is dominated by perturbations of the perimeter of the 
holes, an example of which is shown in Fig. 1(c); this is 
also consistent with previous models used to successfully 
fit experiments''' 9 . We take the radial perturbation Ar 
to be a Gaussian random variable with a mean of zero 
and a standard deviation of a. Two radial perturbations 
are correlated by (Ar(&)Ar($)) = a 2 eT^-VM 1 * 5 itj , 
where the subscript indexes the holes, <pi is the angu- 
lar position of the point measured about the centre of 
the hole, r is the hole radius, and l p is the correla- 
tion length measured around the circumference. The 
change in dielectric constant can be written as Ae(r) = 
(e 2 — £i) Ar 5(^Jxj + yf — r) 6(/i/2 — \z\), where Xi and 
Hi are the components of r measured from the centre of 
the hole, and the Heaviside step function restricts the 
disorder to a slab of thickness h. In this analysis, below, 
we used representative disorder parameters of a — 3 nm 
and l p = 40 nm, which are typical for fabricated samples. 

The failure of the weak contrast model, namely Pdis = 
AeEo is most obvious by considering a small dielectric 
sphere introduced in a homogeneous background with di- 
electric constant £2- It is well known that the weak con- 
trast polarizability Ae must be replaced with the correct 
polarizability 3Ae/(3e2 + Ae) (as can easily be proven 
from Eq. (1) ). Due to local field effects, the macro- 
scopic scattering depends on the microscopic geometry 
of the scatter, and this is also true for disordered waveg- 
uides. This issue has been partly investigated by Johnson 
et al. 15 , who give the correct disorder polarization den- 
sity due to a disorder element at r' as (lu is implicit) 

P(r) = f£l±£ia||E|| (r) + e(r) 7± T> ± (r)j AVS(r - r'), 

(2) 

where ay and 7^ are polarizabilities for the disorder cl- 
ement, e(r) takes a different value depending on which 
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side of the interface r' is located, and AV is the vol- 
ume of the disorder element. Due to the e(r) factor, 
E*(r) ■ P(r) will have terms proportional to |E|| (r)| 2 and 
|D^(r)| 2 and these fields are well defined at a dielectric 
interface. This improved disorder model hides a great 
deal of complexity in the polarizabilities a\\ and 7j_. The 
polarizabilitics arc different for a positive bump (ei ex- 
tending into £2, Q!|| = Qij| ) and a negative bump (a^). 
Further, in general, the polarizabilities are asymmetric 
so that ajj~ ^ — a.7. Finally, the exact polarizabilities 
depend on the precise geometry of the disorder element 
and must be calculated numerically. It is thus no surprise 
that only a single isolated bump was treated in Ref. 13. 

One method to simplify the treatment of the 
polarizabilities is to assume a structure for the 
disorder where the polarizability is known. In 
Ref. 15, a polarization density of the form P(r) — 

Ae(r) ^E|| (r) + e(r) J 5{r— r') was used, namely the 

slowly-varying surface approximation. This takes care of 
the issue of the E(r) being ill-defined at the interface but 
it has some subtle and questionable assumptions. Im- 
plicit in the derivation of this term is the assumption 
that the disordered surface is smooth; effectively that 
the hole radius changes but remains nearly circular. As 
a model for the disorder in real fabricated samples, this 
is rather suspect, and it is known that the interface fluc- 
tuations vary rapidly (cf. Fig. 1(c) and Ref. 20). Our 
solution to this general problem is to use a representa- 
tive model for the polarizability, namely Eq. 2 with the 
numerically calculated polarizabilitics for a cylindrical 
bump 15 , in combination with our numerically generated 
disordered profile. We use this straightforward and qual- 
itative approach to provide essential physics insight into 
optical scattering phenomena that may occur when the 
polarizability is correctly treated, which, as mentioned 
previously, is lacking in all previous disorder models. 

It is easiest to compare these polarization density 
models above using our incoherent scattering calculation 
where the band-edge resonances are suppressed. The 
incoherent average power loss, which is dominated by 
backscatter loss for slow light PC waveguides, is ,! 

(«bac k )= f f drdr' <[E(r)-P(r)] [E*(r')-P*(r')]) > 

(3) 

where the integration is performed over a single unit 
cell and E(r) = ek(r)e lkx . The incoherent averaged 
backscatter loss is plotted in Fig. 2 for the weak contrast 
(blue), smooth surface (green), and cylindrical bump 
(red) models. Although the exact loss predictions differ, 
all the models predict similar loss trends and their mag- 
nitudes are similar. For the roughness statistics typical of 
PC waveguide systems, we find reasonable agreement be- 
tween the backscatter loss predictions of the three mod- 
els, though the smooth surface and weak contrast models 
underestimate the scattering loss. But this by itself, al- 
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FIG. 2: (Color online) Incoherent averaged back scatter loss 
for a single unit cell using the weak contrast (blue), smooth 
surface (green) and cylindrical bump (red) polarization den- 
sity models. For reference, the group index (magenta, dot- 
dash) is shown on the right scale. 

though an improvement, is not a drastic revision of pre- 
vious disorder models. 

As discussed previously, one of the effects neglected 
when one ignores the Go(r,r';u) contribution of the 
Green function is disorder-induced frequency shifts. 
This-hitherto-theoretically unknown phenomenon for 
PC waveguides is closely related to the selection of a suit- 
able polarization model, as the predicted frequency shift 
is sensitive to this term. The first-order mean frequency 
shift due to disorder is given as 21 

(Aa,) = -|/dr<E*(r)-P(r)) (4) 

with E(r) normalized as before. The frequency shift is 
usually taken to be zero for zero-mean surface perturba- 
tions, but as shown previously for a PC cavity 11 ', cor- 
rectly treating local field effects yields a non-zero first 
order frequency shift. Considering the above polariza- 
tion densities, both the weak contrast and smooth sur- 
face models predict (Auj) = due to the symmetry in the 
polarization for positive and negative bumps. The cylin- 
drical bump polarizability predicts a non-zero (Aw), but 
the details differ from Ref. 16 (the simpler PC cavity 
case) due to differing disorder models; for the waveguide, 
we are also dealing with a continuous mode rather than 
a discrete resonance. When using the polarization den- 
sity of Eq. 2, care must be taken when evaluating the 
expectation value since the value of the polarizabilitics 
depends on the direction of the bump. In addition to the 
mean frequency shift, which in general may or may not 
be zero, the RMS frequency shift (Auj)rms = \J (Auj 2 ) 
is certainly not. It is calculated in a similar way to the 
corrected backscatter loss, from 

(Aw 2 ) = ^ J! dr dr' <[E*(r) • P(r)] [E*(r') ■ P(r')]> • 

(5) 

We stress that neither of these two frequency shifts have 
been predicted nor calculated for PC waveguides, yet 
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FIG. 3: (Color online) Mean frequency shift (dashed) and 
RMS frequency shift (solid) using the weak contrast (blue), 
smooth surface (green) and cylindrical bump (red) polariza- 
tion density models. Only the cylindrical bump model has a 
non-zero mean frequency shift. For reference, the group index 
(magenta, dot-dash) is shown on the right scale. 



clearly they are just as important, if not more so, as 
computing the disorder-induced power loss. Figure 3 
plots the mean (dashed) and RMS (solid) frequency shifts 
for the weak contrast (blue), smooth surface (green) 
and cylindrical bump (red) polarization density mod- 
els. The weak contrast and smooth surface models — 
incorrectly predict zero mean frequency shifts and thus 
are not shown. The disorder-induced frequency shifts 
are particularly important for understanding experimen- 
tal transmission spectra. Typically, very near the band 
edge, there will be an abrupt drop in the transmission 
associated with a local frequency shift causing the band 
to be shifted such that the injected frequency is below 
the new band edge. As can be seen from the RMS fre- 
quency shifts, this will be an issue regardless of the po- 
larization model, but can be exasperated by a non-zero 
mean frequency shift. The inset to Fig. 1(a) depicts the 
computed band structure with disorder" - (grey shaded 



band), showing that getting too close to the mode edge 
is actually impossible. Qualitatively, for waveguide mode 
frequencies within a few standard deviations of the mode 
edge (slow light regime) , we expect that somewhere along 
the waveguide, total reflection will inevitably occur. 

These results arc directly relevant to a wide range of 
enhanced light-matter interaction physics that occurs in 
nanophotonic waveguides. Apart from the consistent ob- 
servation of a sudden, dramatic reduction in transmission 
in PC waveguide near the mode edge (underestimated by 
current theories), we cite two other examples: (i) obser- 
vation of the enhanced spontaneous emission with sin- 
gle photon emitters (quantum dots) near the mode edge 
indeed observe a broadened and substantially reduced 
density of states than is expected from the disorder-free 
band structure 23 ; (ii) recent work by Morichetti et al. 
demonstrate that backscattcring is also one of the most 
severe limiting factors in state-of-the art silicon on insula- 
tor nanowires, and our conclusions are directly applicable 
and supportive of these measurements as well: broaden- 
ing and band structure restructuring near the mode edge 
will, inevitably, result in dramatic backscattcring. 

In summary, we have described a theory of disorder- 
induced scattering to include the influence of local field 
effects and high-indcx-contrast perturbations. Our cal- 
culations are shown to increase the predicted waveguide 
losses and result in significant and surprising (asymmet- 
ric) disorder-induced resonance shifts. The band struc- 
ture broadening (schematically shown in Fig. 1(a)) offers 
fresh and important insights into the fabrication limits 
of slow light propagation in PC waveguides and should 
serve as a further warning that propagation modes near 
the mode edge will have an increasingly better chance of 
being completely reflected. Our predictions are consis- 
tent with the wide range of experiments in the literature. 
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